Non-equilibrium emission of complex fragments from p+Au collisions at 2.5 GeV 

proton beam energy 
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Energy and angular dependence of double differential cross sections d 2 cr/dfidE was measured for 
reactions induced by 2.5 GeV protons on Au target with isotopic identification of light products (H, 
He, Li, Be, and B) and with elemental identification of heavier intermediate mass fragments (C, N, 
O, F, Ne, Na, Mg, and Al). It was found that two different reaction mechanisms give comparable 
contributions to the cross sections. The intranuclear cascade of nucleon-nucleon collisions followed 
by evaporation from an equilibrated residuum describes low energy part of the energy distributions 
whereas another reaction mechanism is responsible for high energy part of the spectra of composite 
particles. Phenomenological model description of the differential cross sections by isotropic emission 
from two moving sources led to a very good description of all measured data. Values of the extracted 
parameters of the emitting sources are compatible with the hypothesis claiming that the high energy 
particles emerge from pre-equilibrium processes consisting in a breakup of the target into three 
groups of nucleons; small, fast and hot fireball of ~ 8 nucleons, and two larger, excited prefragments, 
which emits the light charged particles and intermediate mass fragments. The smaller of them 
contains ~ 20 nucleons and moves with velocity larger than the CM velocity of the proton projectile 
and the target. The heavier prefragment behaves similarly as the heavy residuum of the intranuclear 
cascade of nucleon-nucleon collisions. 

PACS numbers: 25.40-h,25.40.Sc,25.40.Ve 

Keywords: Proton induced reactions, spallation, fragmentation 

I. INTRODUCTION of the intranuclear cascade of nucleon-nucleon collisions 

followed by fragment evaporation from excited remnant 
nucleus in competition with fission process. Whereas the 
low energy part of spectra - up to 60 - 80 MeV - seems to 
be reasonably well described by this conventional mecha- 
nism, the high energy part of spectra is generally strongly 
underestimated by any of the existing models. It is more 
flat than the low energy part of the spectrum and its 
slope increases monotonically with the emission angle. 
This behavior indicates the necessity to include non- 
equilibrium processes in the description of the reaction 
mechanism. Authors of Rcfs. [KH12| propose the sur- 
face coalescence of emitted nucleons as process respon- 
sible for high energy part of the 2,3 H and 3 He spectra. 
They claim, however, that such a mechanism is ruled out 
for 4 He and heavier ejectiles [HI EH ■ 

In the present study the task was undertaken to mea- 
sure double differential cross sections {d 2 <j/d£ldE) with 
isotopic identification of the light reaction products from 
proton - gold collisions at proton beam energy of 2.5 GeV, 
extending the range of detected ejectiles to heavier than 
those from previous reported investigations [lol . [Tl| . 

'Electronic address: ufkamys@cyf-kr.edu.pl [correspondin g author] | It should be emphasized that for the gold target the 



The mechanism of proton - nucleus interactions at GeV 
energies is still not well understood. Even the gold nu- 
cleus which is the most frequently studied target, at least 
as concerns the measurements of total cross sections for 
emission of different products (cf. Refs.[ll-Q and refer- 
ences herein) , reveals unexpected phenomena when more 
exclusive observables are investigated. Recently measure- 
ments of differential cross sections in 4-7T geometry were 
undertaken [h], |Tl[ for light charged particles (LCP's), 
i.e. H and He ions, as well as for intermediate mass 
fragments (IMF's) - Li and Be ions. The measurements 
were done at 1.2 GeV proton energy for 1,2 ' 3 H, 3 ' 4,6 Hc, 
6 ' 7 - 8 ' 9 Li and 7 < 9 > 10 Be isotopes and at 2.5 GeV for 

l,2.3 Hj 3,4 H6j and 6,7 Li isotopes [ll| . 

It was found that the shape of energy spectra of emit- 
ted composite particles as well as their angular depen- 
dence cannot be explained using the conventional picture 
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double differential cross sections (d 2 a / (KldE) of interme- 
diate mass fragments, i.e. fragments with mass number 
Ap > 4, were not measured up to now with isotopic 
identification . The only available data are published 
by Letourneau et al. [nj for 6 ' 7 Li. Low statistics of 
isotopically identified data in publication of Herbach et 
al. [H| did not allow to analyze double differential cross 
sections (d 2 a/d^ldE) for isotopes heavier than 4 He. For 
individual isotopes only the analysis of angle integrated 
(da / dE) spectra or energy integrated (da/dQ) angular 
distributions was possible. 

The goal of the present study was to gain new exper- 
imental information on the proton - gold interaction at 
proton energy of 2.5 GeV. Those new double differential 
data (d 2 <r/dildE) should allow to gain deeper insight in 
the mechanism of non-equilibrium processes. 

Details of the experimental procedure are discussed in 
the second section and the obtained data in the next, 
third section. The fourth section is devoted to model 
description of the measured spectra. The interpreta- 
tion, summary and conclusions are presented in the 
last section. Formulae applied in the phenomenological 
parametrization are collected in the appendix. 



II. EXPERIMENTAL PROCEDURE 

The experiment has been performed using the internal 
beam of the Cooler Synchrotron COSY in the Research 
Center in Julich. Due to multiple passing of the circulat- 
ing internal beam through the target it was possible to 
achieve as high luminosity as that which can be reached 
only with very intensive external beam of accelerators 
(with the particle current of order of mA). Circulation 
of the beam without its immediate absorption demanded 
using of very thin, self supporting targets (of order of 
300 fig/ cm 2 ) what in turn resulted in negligible distor- 
tion of the reaction product spectra by interaction of the 
emitted particles with the target. Furthermore, during 
each cycle of injection and acceleration, the protons were 
circulating in the COSY ring slightly below the target, 
being slowly bumped onto the target until the beam was 
completely used up. Then a new cycle was started. The 
speed of the vertical shift of the proton beam was con- 
trolled by feedback of the observed reaction rate to avoid 
overloading of the data acquisition system. 

The scattering chamber and the detecting system was 
described in detail in ref. [3]. There, however, the main 
interest was focused only on performance of the gridded 
ionization chambers which are used in the experiment for 
charge identification as well as for energy measurement 
of the reaction products by means of the Bragg curve 
spectroscopy. For this reason, a description of the other 
components of the detection system, relevant to the data 
discussed in the present paper, will be given here in a 
more detailed way. 

The PISA apparatus consists of nine independent de- 
tection arms comprising various kinds of detectors. Two 



of these arms (placed at 15 and 120° angles in respect to 
the beam direction) are equipped with the Bragg curve 
detectors (BCD), which permit the ^-identification of 
the reaction products and determination of their kinetic 
energies with low detection energy threshold (of about 
1 MeV/nucleon). The telescopes consisted of silicon de- 
tectors are installed at the detection angles of 35, 50, 80 
and 100°. The detectors operate in the ultra high vac- 
uum (UHV) of the COSY accelerator and are cooled to a 
temperature of -10° C. The cooling improves the energy 
resolution of the detectors, thus the isotopes of all ejec- 
tiles up to carbon can be identified. Due to geometrical 
constraints the silicon telescope detectors placed in the 
vacuum at 35, 50, and 80° cannot be supplied with addi- 
tional detectors. Consequently light charged particles of 
high energies, not being stopped in the silicon detectors, 
cannot be detected. The upper limit of energies of par- 
ticles stopped by these telescopes is about 30 - 40 MeV 
for protons, deuterons and tritons but increases signif- 
icantly for heavier particles, e.g. for alpha particles it 
is around 120 MeV. Therefore these telescopes are suit- 
able to measure the low energy part of the spectra for 
hydrogen isotopes, a large range of energies for helium 
isotopes and the full energy spectra of intermediate mass 
fragments. The silicon detector telescope placed at 100° 
in the ultra high vacuum has another construction than 
the telescopes mentioned above, thus it was possible to 
supplement it with a 7.5 cm thick Csl scintillator detec- 
tor standing behind it in the air (outside the chamber), 
separated by a steel window of 50 /im thickness from the 
ultra high vacuum of COSY. At three angles (15.6, 20, 
and 65°), the telescopes built of silicon detectors with 
7.5 cm Csl scintillator detectors standing behind them 
are positioned outside the chamber. The destination of 
these telescopes as well as that at 100° is to significantly 
increase range of energies of detected light charged par- 
ticles and IMF's. 

The particles observed at different angles in the present 
experiment and the energy ranges covered by the detect- 
ing system are listed in the Table [J and in the Table ITT1 for 
isotopically identified and elementally identified cjectilcs, 
respectively. 

The absolute normalization of the data was achieved 
by comparison of the total cross section for 7 Be ejec- 
tiles extracted from angular and energy integration of 
the spectra measured in the present experiment with the 
cross section obtained from parametrization of experi- 
mental 7 Be production cross sections in proton-nucleus 
collisions, ref. Accuracy of the absolute normaliza- 

tion was estimated to be better than 10%. 



III. EXPERIMENTAL RESULTS 

In the present study the double differential spectra 
d 2 a/dftdE were determined for the first time for many 
isotopically identified intermediate mass fragments emit- 
ted from proton - gold collisions at GeV energies. This 
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TABLE I: Range of energies (in MeV) of isotopically identified reaction products detected at various scattering angles 
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TABLE II: Range of energies (in MeV) of elementally identified reaction products detected at various scattering angles 
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50.5 - 116.5 


48.5 - 102.5 


12.5 - 57.5 


N 


14.5 - 69.5 


56.5 - 116.5 


48.5 - 108.5 
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28.5 - 98.5 


Al 
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concerns 6 He, 8,9 Li, 7 ' 9,10 Be, 10 ' n ' 12 B spectra which were 
not measured by Letourneau et al. at 2.5 GeV [TlT | 
whereas the experiment of Herbach et al. [l(| at 1.2 GeV 
which detected IMF's lighter than boron had statistics 
allowing to extract only elemental spectra. Typical spec- 
tra for isotopically identified particles from the present 
experiment are shown in Fi g. [1] together with data mea- 
sured by Letourneau et al. [llj. Excellent agreement of 
the present data with those published by Letourneau et 
al. was achieved for all products measured in both exper- 
iments, i.e. 1,2,3 H, 3,4 He, and 6,7 Li. Note, that statistical 
errors, which are only shown for selected 7 Li data of Ref. 
[ill ], present indeed typical errors for all 7 Li data from 
that paper. 

The energy distributions of emitted ejectiles have 
shapes resembling Maxwcllian evaporation spectra, but 
because of an instrumental low-energy cutoff it was not 
possible to observe the maxima of these distributions for 
fragments heavier than boron. Although for heavier frag- 
ments the Coulomb barrier moves the position of maxi- 
mum of the yield towards higher values, the large energy 
loss in first silicon detector of telescope prevent us from 
detecting heavy ejectiles at energies close to the maxi- 
mum of the energy distributions. To avoid this problem, 



i.e. to measure low energy part of the spectra, two Bragg 
curve ionization chambers (BCD) were applied. They 
were placed at 15° and 120° scattering angles. Since 
BCD's allow for the identification of the charge of ejec- 
tiles only, the spectra of heavy products, i.e. C, O, N, 
F, Ne, Na, Mg and Al were obtained only with elemental 
identification. As far as we know, similar spectra were up 
to now investigated for the gold target only in the exper- 
iment performed at 1.0 GeV by Kotov et al. [It}, where 
double differential cross sections were measured without 
isotopic identification. 

Typical properties of the spectra, characteristic for 
all ejectiles, are depicted in Fig. [TJ At low energy the 
angle independent - Maxwcllian like - contribution is 
well visible. This isotropic energy distribution may be 
reproduced by the two stage model discussed below. The 
dashed line shown in the figure, represents predictions 
of this model. Another contribution, i.e. an exponential 
distribution, strongly varying with angle is present at 
higher energy in all experimental spectra. The slope of 
this anisotropic energy contribution increases with the 
angle, what may be interpreted as effect of fast motion 
of an emitting source in the forward direction. 
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FIG. 2: Open circles represent experimental energy spectrum 
of 7 Li particles measured at 35° in the current experiment. 
The lines present result of phenomenological model described 
below; the short-dashed line shows contribution of the slow 
emitting source, the long-dashed line depicts contribution of 
the fast source, whereas the solid line presents sum of both 
contributions. Please note, that the shape of this experimen- 
tal energy distribution as compared with spectra shown in 
right part of Fig. [1] also confirms the monotonic change of 
the exponential slope with the scattering angle. 



FIG. 1: Typical energy spectra of 4 He (left column) and 7 Li 
particles (right column) measured in the present experiment 
(open circles) and published in Ref. [tP ] (full dots) for corre- 
sponding emission angles. The lines show prediction of evap- 
oration of 4 He and 7 Li evaluated by means of GEM program 
of S. Furihata [l5|, [TH from excited residual nuclei of the first 
stage of the reaction with properties extracted from BUU cal- 
culations. 



IV. THEORETICAL ANALYSIS 

In the most frequently considered scenario of the 
proton-nucleus collision at GeV proton energies it is as- 
sumed that reaction proceeds via two stages. 

In the first stage of the reaction the proton impinging 
on to the target nucleus initiates a cascade of nucleon- 
nucleon collisions which leads to emission of several fast 
nucleons and pions, and to excitation of the nucleus. 
This fast stage of reaction is described by intranuclear 
cascade (INC) model, e.g. [HJll, Boltzmann-Uchling- 
Uhlcnbcck (BUU) model, e.g. [2fJ] or by quantum molec- 
ular dynamics (QMD) model, e.g. j21j. The first of 
the mentioned models gives an account of the nucleon- 
nucleus interaction by static (time-independent) mean 
field, the BUU allows for dynamic evolution of the mean 
field as caused by time dependence of an average nucleon 
density, and the QMD treats the nucleon-nucleus interac- 
tion as a time dependent sum of elementary two-nucleon 
and three-nucleon interactions of all nucleons. The QMD 



introduces the largest fluctuations of the density distri- 
bution of nucleons and, therefore, allows for emission of 
clusters of nucleons from the first stage of the reaction. 
The static mean field description used by INC model au- 
tomatically precludes possibility of nucleon distribution 
fluctuations. The BUU model takes into consideration a 
time dependent modification of the nucleon density dis- 
tribution, however, the averaging over many test parti- 
cles, present inherently in BUU, prohibits appearing of 
fluctuations large enough for nucleon clusters emission. 
The emission of fast nucleons (in the case of INC and 
BUU) or fast nucleons and light clusters (in the case of 
QMD) terminates after a short time, leaving the residual 
excited nucleus in a status close to the thermodynamic 
equilibrium. 

The second stage of reaction consists in the evapo- 
ration of nucleons and clusters from this equilibrated 
system, which can also undergo fission with emission 
of two heavy fragments. Thus, in the two-step model 
of reaction mechanism, the non equilibrium emission of 
nucleons and clusters can appear only in the first stage 
of the reaction. It is believed that statistical model codes 
like, e.g. GEM 03 or GEMINI [H are capable to 
well reproduce emission of nucleons and fragments from 
equilibrated, excited nucleus. Therefore, observation of 
any disagreement of the data with predictions of the 
two-step model would suggest that (i) the model is 
not adequate to the real situation (e.g. an additional, 
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FIG. 3: Open circles represent typical spectra of protons, 
deuterons and tritons measured in the present experiment 
by telescope consisted of silicon semiconductor detectors and 
7.5 cm thick scintillating detector Csl placed at scattering 
angle of 65 degree in respect to the proton beam. The dashed 
lines show evaporation contribution evaluated by means of 
the BUU and Generalized Evaporation Model whereas the full 
lines correspond to phenomenological model of two emitting 
sources described below. Note change of the scale for the 
triton spectrum. 



1 




To 10 



CM 

T3 




0,01, ' 



1E-3I I 

30 60 90 120 150 

E [ MeV ] 



FIG. 4: Typical energy spectra of helium ions 3,4,6 He mea- 
sured in the present experiment by telescope consisted of sil- 
icon semiconductor detectors placed at scattering angle of 35 
degree in respect to the proton beam - open circles. Note 
different vertical scales for each spectrum. The lines have the 
same meaning as in Fig. 



A. Boltzmann-Uehling-Uhlenbeck model and 
evaporation model 



intermediate stage of the process is necessary before 
achieving thermodynamic equilibrium), or (ii) descrip- 
tion of the emission of particles from the first stage of 
the reaction (nucleons or clusters) is not properly taken 
into consideration. 



The present data were compared with results of a two 
stage model in which the Boltzmann-Uehling-Uhlenbeck 
transport equation [20| has been applied for the descrip- 
tion of the first step of the proton - nucleus collision 
leading to emission of fast nucleons leaving the heavy 
excited remnant in a state close to equilibrium. Monte 
Carlo computer program developed at Giessen Univer- 
sity [IK was utilized to simulate this first stage of the 
reaction and to find properties of excited residual nuclei. 
Deexcitation of these nuclei, which proceeds by emission 
of nucleons and complex fragments, was calculated in 
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the frame of statistical model using the GEM (General- 
ized Evaporation Model) computer program of Furihata 
flU Theoretical energy spectra of various ejectiles 
found from this two stage model are shown in Figs. [T]-[5] 
as dashed lines. It can be concluded from examination of 
these figures that the model predictions describe well low 
energy part of spectra for hydrogen, helium and lithium 
isotopes. For heavier ejectiles the theoretical cross sec- 
tions underestimate the experimental data. Moreover, it 
can be observed that the high energy part of the spectra 
is clearly not reproduced by the two stage model, which 
predicts much steeper slope of the spectra than is ob- 
served experimentally. Thus, another mechanism seems 
to give a significant contribution to the proton - nucleus 
reactions. As concerns hydro gen and helium production, 
the authors of [lj| , [H| , and [l2| papers propose the co- 
alescence of nucleons as the mechanism responsible for 
this effect, however, no microscopic model is able to re- 
produce observed effects for heavier composite ejectiles. 

An extensive comparison of predictions resulting from 
the models mentioned above with our experimental data 
presented here will be subject of a forthcoming paper. We 
restrict ourselves here on conclusions we can draw from 
the application of a phcnomcnological model described 
in the next section. 

Following properties of the spectra should be taken 
into consideration when looking for an appropriate phe- 
nomenological model: 

(i) The position of the peak present at low energies in 
the experimental spectra of all observed particles 
(and its height for light ejectiles) is quite well re- 
produced by the two stage model discussed above. 
This means that the mechanism described by this 
model gives a large contribution to the reaction and 
therefore it must be taken into account in the frame 
of any phenomenological model. 

(ii) The slope of the exponential, high energy tail of the 
experimental spectra for all ejectiles varies mono- 
tonically, increasing strongly with the scattering 
angle as can be seen from Fig. [1] Such a behav- 
ior is in contradistinction to properties of the spec- 
tra evaluated in the frame of the two stage model, 
which are almost independent of angle. This in- 
dicates that high energy particles are not emitted 
from heavy residuum of the intranuclear cascade 
but from another source which moves much faster 
than the residuum. 

These arguments call for using of a phenomenological 
model of two emitting sources; one source moving slowly 
would imitate emission from heavy residuum of the in- 
tranuclear cascade whereas the second source should sim- 
ulate emission from faster (and thus probably lighter) 
nuclear system. Of course, one could imagine that more 
than two sources of emitted particles are necessary for 
reasonable description of the data. The applied model of 
two moving sources corresponds to minimal number of 
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FIG. 5: Typical spectra of lithium ions 6,7 ' 8 ' 9 Li measured in 
the present experiment by telescope consisted of silicon semi- 
conductor detectors placed at scattering angle of 35 degree 
in respect to the proton beam - open circles. Note different 
vertical scales for 6,7 Li and 8,9 Li. The lines have the same 
meaning as in Fig. [3] 



parameters necessary to fulfill qualitative demands put 
on the model by the experimental data. 



B. Phenomenological model of two moving sources 

In the frame of the phcnomcnological model of two 
moving sources the angular and energy dependence of the 
double differential cross sections d 2 a/dVldE is described 
by analytical formulae. The details of the model and 
interpretation of its parameters are presented in the Ap- 
pendix. An example of the description of the experimen- 
tal energy spectrum by the two source model is shown 
in Fig. [5J The symbols depict the data from present ex- 
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FIG. 6: Typical spectra of beryllium ions 7,9,10 Be measured in 
the present experiment by telescope consisted of silicon semi- 
conductor detectors placed at scattering angle of 35 degree in 
respect to the proton beam - open circles. The lines have the 
same meaning as in Fig. [3] 



FIG. 7: Typical spectra of boron ions 10 ' n ' 12 B measured in 
the present experiment by telescope consisted of silicon semi- 
conductor detectors placed at scattering angle of 35 degree in 
respect to the proton beam - open circles. The lines have the 
same meaning as in Fig. [3] 



periment obtained for 7 Li ejectiles detected at scattering 
angle 35° whereas the lines show result of the fit of the 
phenomenological model. The short-dashed line presents 
contribution from the slowly moving source, the long- 
dashed line shows contribution from the fast source and 
the solid line corresponds to sum of both contributions. 
As can be seen, very good description of full energy spec- 
trum could be achieved. 

The parameters of the theoretical formula of the two 
moving sources model have been searched for by fitting 
simultaneously experimental spectra at several scattering 
angles for each ejectile. Exceptions from this rule were 
the spectra of ejectiles heavier than oxygen (F, Ne, Na, 
Mg. and Al). which were measured only at these two an- 
gles at which Bragg curve ionization chambers have been 



positioned, i.e. at 15° and 120°. Such spectra were fitted 
assuming that only one moving source gives contribution 
to the reaction. Furthermore, the spectra of C, N, and 
O which were measured both, by silicon detectors at 35°, 
50°, 80°, and 100° as well as by Bragg curve ionization 
chambers at 15° and 120° degree were fitted using one 
emitting source and two emitting sources. The parame- 
ters of sources for light charged particles and isotopically 
identified IMF's are listed in Table Hill whereas those for 
heavier IMF's, which were only elementally identified, are 
collected in Table JV] 

The first source should simulate evaporation of par- 
ticles from heavy remnant of the first stage of the re- 
action, i.e. intranuclear cascade of nucleon-nucleon col- 
lisions. Thus, its velocity was fixed at value ,3=0.002 
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FIG. 8: Typical spectra of carbon, nitrogen and oxygen ions 
measured in the present experiment without isotopic separa- 
tion by telescope consisted of silicon semiconductor detectors 
placed at scattering angle of 35 degree in respect to the pro- 
ton beam - open circles. The lines have the same meaning as 
in Fig. [3] 



(in units of velocity of light) as it was extracted from 
BUU calculations. This value was constant for all cal- 
culations. Other parameters characterizing the source, 
i.e. fc-paramctcr (ratio of the actual height of Coulomb 
barrier to its value found from simple estimation for two 
touching, charged spheres), T-parameter (apparent tem- 
perature), and a (energy and angle integrated cross sec- 
tion for production of given ejectile) were free parameters 
of the fit. 

All parameters of the second source were freely mod- 
ified in fits since no hypothesis concerning origin of this 
source was made before the analysis. Usually the pro- 
gram was able to find unambiguously the best parame- 
ters, corresponding to the minimum value of chi-squarc. 
In such a situation the routine provides estimation of er- 
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FIG. 9: Examples of energy spectra for heavier, elementally 
identified IMF's obtained in the current experiment by means 
of Bragg curve ionization chambers. The open circles repre- 
sent experimental data and the solid lines show results of 
phenomenological model analysis obtained with assumption 
of only one moving source, discussed in the next section of 
the paper. 



rors. In some cases, however, the valley of chi-squarc 
values was so complicated that the program was not able 
to produce reasonable estimation of errors. The ambigu- 
ity of parameters lead sometimes the searching procedure 
to nonphysical values of the parameters, as e.g. negative 
height of the Coulomb barrier. Then it was necessary to 
fix these parameters at values, which still have physical 
meaning. Such values of parameters are quoted in the 
tables as closed in square parentheses. 

Thorough inspection of the parameter dependence 
leads to the following conclusions: 

1 . The contribution o\ of the first ( slow ) emit- 
ting source is comparable to contribution 02 
of the second source. 

It is illustrated by Fig. [TQ] in which ratio of the 
total cross section for emission of ejectiles from the 
first source to the sum of the total cross sections for 
emission from both sources is shown as function of 
mass number of the ejectiles. The average value of 
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TABLE III: Parameters of two moving sources for isotopically identified products 
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U.U-tU ZLZ U.UU/ 


91 fi 4 1 4 
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1.0 It 1.1 


9 1 
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U.OU ZLZ U.U'tfc 


H I _j- 0.8 


9^ 14-17 
zo.o m -L.i 


n 44 4 n ri4 

U.44 It U.Ul 


0^4 4- 00^ 

U.UOI ZLZ U.UUO 


91 7 4 fi 

ZO. i It u.u 


14 K 4- 1 7 
11.0 it 1.1 


9 
z.u 


Y Ti 


fiQ 4- n 0^ 

U.OO ZLZ U.UO 


11.U it u.u 


c:n 8 4 9 fi 
ou.o m z.u 


(i ¥ 4 n ni 

U.OU It U.UO 


oq^ -1- 00^ 

U.UOO ZLZ U.UUO 


on Q 4 f) ^ 
zu.y it u.o 


on 1 4- 9 fi 
ZU.O It z.u 


0.1 




Q0 + 0Q 


11.9 ± 1.5 


9.1 ± 1.4 


45 4 05 


09Q 4- 00^1 

U.UZiC/ —1— U.UUO 


18.0 it 1.0 


6.4 ib 1.5 


9 1 
z . 1 


u Li 


1.00 ± 0.22 


10.4 ± 3.0 


2.1 ± 0.5 


0.39 ± 0.07 


0.025 ± 0.003 


18.2 ib 1.6 


2.1 ib 0.6 


1.2 


'Be 


0.92 ± 0.27 


11.2 ± 4.3 


2.6 ± 0.8 


0.48 ± 0.05 


0.038 ± 0.005 


24.0 ± 1.2 


4.6 =b 0.9 


1.4 


y Be 


0.86 ± 0.12 


9.6 ± 1.7 


12.5 ± 1.9 


0.53 ± 0.06 


0.020 ± 0.005 


16.6 =b 0.8 


8.1 ± 2.3 


1.4 


lu Be 


0.90 ± 0.08 


11.8 ± 1.2 


10.0 ± 1.4 


0.44 ± 0.04 


0.026 ± 0.004 


14.5 ± 0.9 


6.8 =b 1.5 


1.3 


iu B 


0.85 ± 0.20 


10.5 ± 3.4 


6.6 ± 1.3 


0.73 ± 0.14 


0.020 ± 0.010 


18.2 ± 2.7 


2.7 ± 1.7 


1.8 


n B 


0.93 ± 0.18 


10.5 ± 2.1 


12.8 ± 2.5 


0.50 ± 0.05 


0.022 ± 0.004 


14.5 =b 0.7 


12.8 ± 2.8 


1.7 




0.87 


8.8 


1.6 


0.73 


0.012 


13.2 


5.1 


1.0 



TABLE IV: Parameters of one or two moving sources for elementally identified products 





Slow source 


Fast source 




Ejectile 


ki 


Ti/MeV 


cTi /mb 




02 


r 2 /MeV 


a-2/mb 


x 2 


C 


0.879 


12.3 


28.4 


0.150 


0.0367 


15.8 


11.8 


3.34 










0.759 =b 0.032 


0.0076 =b 0.0007 


13.4 ± 0.5 


33.4 ± 0.9 


3.26 


N 


1.00 


12.2 


15.9 


0.206 


0.0382 


14.5 


9.6 


1.47 










0.822 ± 0.045 


0.0118 =b 0.0007 


13.4 ± 0.6 


18.6 ± 0.7 


1.41 


O 


0.75 


14.1 


1.3 


0.71 


0.0129 


12.2 


11.8 


0.78 










0.784 ± 0.045 


0.0118 =b 0.0007 


13.5 =b 0.6 


13.1 ± 0.6 


0.75 


F 








0.545 =b 0.079 


0.0105 =b 0.0014 


17.0 ± 3.2 


5.12 ± 0.32 


0.40 


Ne 








0.670 =b 0.098 


0.0090 =b 0.0018 


15.3 =b 1.9 


5.25 ± 0.44 


0.54 


Na 








0.801 =b 0.117 


0.0104 ± 0.0020 


12.0 ± 0.7 


5.99 ± 0.63 


0.61 


Mg 








0.577 =b 0.090 


0.0103 =b 0.0019 


12.0 ± 1.2 


4.59 ± 0.52 


0.59 


Al 








0.772 ± 0.180 


0.0036 =b 0.0013 


10.0 =b 0.7 


5.24 ± 0.91 


0.23 



the ratio a\/(<Ji +a 2 ) is equal to 0.560 ± 0.044. It 
should be, however, emphasized that rather large 
deviations from the average value appear for indi- 
vidual ejectiles. For example, almost 90 % of alpha 
particles is emitted from the slow source, whereas 
it is the case only for ~ 25 % of protons. 

2. The parameters of the slow source have 
values which agree with the assumption 
that this source simulates evaporation from 
a heavy nucleus corresponding, e.g. to the 
residuum of the target after the intranuclear cas- 
cade of nucleon-nuclcon collisions, namely: 

• The apparent temperature of the slow source 
is independent of the mass of emitted inter- 
mediate mass fragments, what can be seen in 
Fig. ITTI where the slope parameter of the solid 
line is equal to zero within the limits of er- 
rors : - 0.15 ± 0.17. Its stability indicates that 
the recoil effect of the source during emission 
of fragments is negligible, thus, mass of the 



source must be much larger than masses of ob- 
served IMF's (cf. Appendix). Moreover, the 
horizontal line: T = 11.9 ± 1.5 MeV, fitted to 
temperature values extracted from spectra of 
IMF's reproduces also quite well values of the 
apparent temperature for light charged parti- 
cles (H and He ions) as can be expected for 
emission from heavy residuum of the intranu- 
clear cascade of nucleon-nuclcon collisions. 

• The k parameter, which determines the height 
of the effective Coulomb barrier between the 
emitted fragment and the rest of the emitting 
source (cf. Appendix) is very close to unity, 
what means that the charge of the source does 
not differ significantly from the charge of the 
target. It is illustrated by Fig. [T2] where k 
parameters for both sources are shown for in- 
dividual ejectiles. The full squares represent 
the slow source and open squares correspond 
to the second, fast source. The dashed line 
k = 1 is shown in the figure to facilitate judg- 
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FIG. 10: Ratio of total production cross section corresponding 
to emission from the slow source to sum of the total cross 
sections representing emission from both sources versus mass 
number of the detected reaction product. The parameters 
<7i and o 2 were taken as total cross sections because they 
correspond to angle and energy integrated double differential 
cross sections d 2 a/d£ldE. The solid line shows average value 
of the ratio (0.560 ± 0.044). 



merit on the magnitude of the k - parameter. 

3. The second (fast) source is much lighter 
than the residual nucleus of the intranuclear 
cascade because: 

• Its velocity is always larger than limiting ve- 
locity of the proton-target center of mass 
which would be obtained only when total 
beam momentum is transferred to the target 
(J3 w 0.018). Fig. [H illustrates this fact show- 
ing /? 2 values for individual ejectiles as well as 
the horizontal line (3 = 0.018. 

• The Coulomb barrier between the source and 
the ejectile is several times smaller than the 
Coulomb barrier of two touching charged 
spheres representing the target nucleus and 
the ejectile. This is well illustrated by open 
squares in Fig. [T2] 

• The recoil effect (cf. Appendix) is clearly visi- 
ble in the dependence of the apparent temper- 
ature of the source on the mass of the ejectile 
as it is shown in Fig. [JTJ- open squares. 

4. The fast source describing LCP's emission 
(hydrogen and helium ions up to A He) is 
much lighter than the fast source responsi- 
ble for emission of intermediate mass frag- 
ments. 



FIG. 11: Apparent temperature of the slow source (full 
squares and solid lines) and that of the fast source (open 
squares and dashed lines) versus mass number of the ejectiles. 
The lines were fitted separately for light charged particles ( 1 H 
- 4 He) and intermediate mass fragments (Aj? > 6). 



This may be inferred from different recoil effects 
visible as different slopes of two lines which de- 
scribe the dependence of the apparent temperature 
on the mass of ejectiles (cf. Fig. ITT]) . The line cor- 
responding to LCP's is more steep, giving the mass 
of the source equal to As~(8 ± 2) nucleons, and 
very high temperature of the source r=(62 ± 7) 
MeV (cf. Appendix for meaning of r). Velocity of 
this light source is very high; (3 = 0.05 — 0.15. 
Such a source can be, perhaps, identified with the 
fireball created by the proton impinging on to the 
target together with nucleons present on its straight 
line way through the target nucleus [24[ , 

The line describing temperature of IMF's corre- 
sponds to mass of the source As=(20 ± 3) nucleons 
and its temperature r= (33 ± 2) MeV. Velocity of 
this source is much smaller (/? = 0.02 — 0.04) than 
velocity of source emitting LCP's. 

Very different properties of the fast source emitting 
light charged particles (LCP's) and the fast source emit- 
ting intermediate mass fragments (IMF's) leads to con- 
clusion that the picture of two sources is oversimplified. 
The presence of a fireball, which can give contribution 
to emission of LCP's only and occurrence of the light 
(As ~ 20), fast source emitting LCP's as well as IMF's 
may be interpreted as indication of a three body decay 
of the target nucleus. The third partner of such a decay 
would be heavy and hardly distinguishable from heavy 
residuum of the intranuclear cascade, therefore its oc- 
currence could be described effectively by the same slow 
source. 
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FIG. 12: The ejectile mass number dependence of the factor 
scaling the Coulomb barrier of two touching spheres to the 
actual height - necessary for good description of the data. 
Full squares represent the slow source and open squares show 
results for the fast source. The line k=l is also depicted to 
facilitate interpretation of the figure. 



FIG. 13: Velocity of the second (fast) source as a function 
of emitted fragment mass number - open squares. Thin solid 
line /?2 =0.018 represents velocity of the common center of 
mass of the proton projectile and the target and thick solid 
line /3i=0.002 shows velocity of the first (slow) source fixed 
at velocity of heavy residuum of intranuclear cascade. 



V. SUMMARY AND CONCLUSIONS 

The double differential cross sections (d 2 a j 'dD.dE) 
were for the first time measured with good statistics 
for isotopically identified intermediate mass fragments 
produced by interaction of 2.5-GeV protons with the 
gold target. The following individual isotopes of the 
elements from hydrogen to boron were resolved: 1,2 ' 3 H, 
3 ' 4 - 6 He, 6 ' 7 < 8 < 9 Li, 7 > 9 ' 10 Be, w > n > 12 B, whereas for heavier 
ejectiles (from carbon to aluminium) only elemental 
identification was done. The energy spectra for all 
nuclear fragments, determined at several scattering 
angles, appear to be of the Maxwellian shape with 
exponential, high energy tail. The low energy part of 
the distribution is almost independent of angle, but 
the slope of high energy tail of the spectrum increases 
monotonically with the angle. The shape of the angle 
independent part of spectra can be reproduced by the 
two-stage model of the reaction, i.e. intranuclear cascade 
of the nucleon-nucleon and meson-nucleon collisions 
followed by statistical emission from an equilibrated 
residual nucleus. However, the absolute magnitude 
of the spectra predicted by two-stage model, using 
Boltzmann-Uehling-Uhlenbeck program [23| for the in- 
tranuclear cascade and Gencralizcd-Evaporation-Modcl 
(GEM) 0, m for statistical emission of fragments, is 
in agreement with the experimental data only for the 
light charged particles (H and He ions). Furthermore, 
the theoretical cross sections underestimate significantly 
the yield of heavier fragments at high kinetic energies 
for all ejectiles. This indicates that another mechanism 



plays an important role besides the standard two-stage 
mechanism. 

To get information on possible origin of this additional 
mechanism a phenomcnological analysis was performed 
assuming that the ejectiles originate from two moving 
sources. The slow moving source was identified with 
the heavy remnant nucleus of the first stage of the 
two-step process mentioned above while no assumptions 
have been made as concerns the second emitting source. 
The properties of both sources were treated as free 
parameters with exception of the velocity of the slow 
source which was taken to be equal to the velocity of the 
heavy residual nucleus from the intranuclear cascade, 
namely ft =0.002. 

Excellent agreement of the phenomcnological 
paramctrization with experimental data was achieved 
with values of the parameters varying smoothly from 
ejectile to ejectile. Their behavior indicates that the 
parameters of the slow source are compatible with 
the assumption that it is a heavy nucleus which may 
be described as equilibrated system of ~ 12 MeV 
apparent temperature, whereas the second source has 
completely different properties. It corresponds cither 
to very small (As ~ 8), very hot (r ~ 62 MeV) and 
fast (0 = 0.05 - 0.15) fireball or to heavier (A s ~ 20), 
colder (t ~ 33 MeV), and slower (j3 = 0.02-0.04) cluster. 

The described above properties of two emitting 
sources observed in the interaction of 2.5 GeV protons 
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with the gold target lead to a conclusion that two 
different mechanisms are observed giving almost the 
same contribution to the cross sections. First of them is 
compatible with the standard, two-stage model whereas 
another one seems to be similar to the picture of cold 
break-up proposed by J. Aichelin, J.Hiifner and R. 
Ibarra 12511. 



For example, it can be expected that similar phenomena 
have to appear for other heavy target nuclei if they are 
observed for the gold target. On the contrary, it is not 
obvious whether proton induced reactions on light tar- 
gets or at significantly different proton energies should 
show similar behavior. 



In the model of cold break-up the energetic proton 
bombarding the target drills a cylindrical hole through 
the nucleus causing that the deformed remnant of the 
collision breaks up into two pieces. Thus, three corre- 
lated groups of nucleons appear after first, short stage 
of the reaction : (i) the fast, small cluster consisted of 
the nucleons which were placed within the cylinder with 
the axis along to the projectile path,(ii) two clusters 
- products of the break up. All three clusters act as 
sources emitting light charged particles, whereas two 
heavier clusters are also responsible for emission of 
intermediate mass fragments. The two latter clusters 
are produced in result of a dynamical process in which 
the "wounded" nucleus cannot come to its ground state 
by emission of nucleons or small clusters, however, the 
correlation of the fast group of nucleons knocked out by 
the projectile is of another, kinematic origin. The high 
energy proton impinging on to the nucleus sees it - due to 
Lorentz contraction - as a narrow disc. Therefore all the 
nucleons which lie on the path of the projectile interact 
simultaneously, as one entity, with the projectile. It was 
shown that this collectivity affects the multi particle 
production in proton-nucleus collisions [26| as well as 
manifests itself in enhanced dependence of momentum 
transfer on projectile energy in deep-spallation reactions 
[27[ | . Thus, it is not surprising that such correlated group 
of nucleons can appear as a hot, fast moving source 
emitting light charged particles observed in the present 
experiment. Of course, it cannot give contribution to 
emission of intermediate mass fragments because the 
source is consisted of a few nucleons only. 

It might be something confusing why in the present 
parametrization only two sources were necessary for 
good description of the data if the postulated cold 
break-up mechanism of the reaction calls for presence of 
three sources. This apparent inconsistency is easy to be 
removed: The slow, heavy source represents the heavy 
residual nucleus produced by the standard two-step 
model of the reaction and/or the heavy fragment from 
the break up . The light, fast source is responsible 
for simulation of the hot fireball (for light charged 
particles) or the lighter fragment from the break up (for 
intermediate mass fragments). 

Consistency of the present phenomenological descrip- 
tion with the cold break up picture of the reaction cannot 
be assumed as a proof of the underlying reaction mecha- 
nism. Additional experimental facts should be searched 
for establishing the confidence in such an interpretation. 



APPENDIX: PHENOMENOLOGICAL 
PARAMETRIZATION 

In this Appendix assumptions and details of the 
formulation of two moving source model are discussed. 
The content of the appendix is very close to information 
contained in the original paper of Westfall et al. [28| . 
however, the additional modification and properties 
introduced in the model need to be discussed for proper 
understanding of the performed analysis. 

The model assumes that the nucleons and composite 
particles are emitted from two moving sources with the 
following properties: 

(i) Each source moves along the proton beam direction, 

(ii) Angular distribution of emitted particles is isotropic 

in the source rest frame, 

(iii) Distribution of the kinetic energy E* available in the 
two-body break up of the source has a Maxwellian 
shape characterized by the temperature parameter 



E* 

T 



dE*dn* 



2(ttt) 



3/2 



The distribution is normalized in such a way that 
integration over angles and energies gives the total cross 
section equal to the parameter a. 

Since the mass of the source As is finite, the energy 
and momentum conservation laws cause that the energy 
E of the observed particle of mass Ap differs from the 
full kinetic energy E* available in the source frame: 



E* 



vE where v 



A, 



A s - At 



thus the energy distribution of the emitted fragment in 
the rest frame of the source is given by: 



dE'dU 



2(ttt) 



3/2 



vE 

T 



This formula is usually applied without explicit writing 
the recoil correction, i.e. by introducing so called appar- 
ent temperature T = t jv: 



d 2 o 
dE'dVL' 



2{-kT) 



3/2 



'W exp 



E 
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Such a form of this formula is used also in the present 
paper. 

It is worth to note, that the recoil of the source gives a 
possibility to extract the information on the source tem- 
perature parameter r as well as on the mass of the source 
As from linear dependence of the apparent temperature 
T on the fragment mass Ap: 

The charged particles emitted from the source 
must overcome the Coulomb barrier, what significantly 
changes the shape of the low energy part of their spec- 
trum. The presence of the barrier may be taken into ac- 
count by shifting the argument in the Maxwell formula 
by the height of the barrier, as it was originally proposed 
in Ref. [28[ or by multiplying the Maxwell distribution 
by the transmission factor. The first method is equiva- 
lent to the application of a sharp cut-off what is a too 
crude approximation, thus the result must be averaged 
over some distribution of heights of the barriers [28| . The 
second method explicitly introduces a smooth variation 
of the transition probability with energy, however, this 
method also must introduce some assumptions concern- 
ing height and curvature of the barrier. In the present 
work, the probability P to overcome the Coulomb barrier 
was parameterized in the following form: 



P = 



1 



1 + exp 



' E-k-B 



)] 



where B is the Coulomb barrier of two touching spheres 
corresponding to the emitted fragment of mass number 
Ap and charge number Zp and to the remaining part 
of the source with the mass number of (As — Ap) and 
charge number (Z$ — Zp): 



D 



Zp(Z s - Z F )e 2 



1.44 (A F 1/3 + (A s - A F ) l /^j ' 



The quantities k and d are the parameters. The first 
parameter (k) gives magnitude of the Coulomb barrier 
in units of B. To avoid ambiguity of B determination 
arising from the fact that at least two different sources 
are present in the current analysis, we evaluated B value 
assuming that Zs = 79 and As = 197, i.e. there are 
atomic and mass numbers of the target. Such value of B 
is a good approximation of the Coulomb barrier for heavy 
residua of the intranuclear cascade, thus one should ex- 
pect that then the k parameter is close to unity. How- 
ever, with such definition of B, the k parameter should 
be significantly smaller than unity for light sources. The 
parameter k was searched by looking for the best fit of 
model spectra to the experimental data. The second pa- 
rameter (d) was arbitrarily fixed in the present analysis 



by keeping constant the ratio of the height of the barrier 
kB to its diffuseness parameter d: kB/d =5.5. 

The explicit introduction of the barrier penetration 
factor P gives finally the following formula for the double 
differential cross section d 2 a/dE dfl : 



exp 



dE'dQ! 
I(B,d,T) 



4wTV*I(kB,d,T) 1 + exp 

oo 

dx ■ \fx ■ exp (—a;) 



kB-E' 



1 + exp ( 



kB — Tx i 



The integral I(B,d, T) used for normalization of the dis- 
tribution (preserving previous interpretation of a pa- 
rameter) has been evaluated numerically by the Gauss- 
Lagucrre method. 

It is necessary to transform the model double differ- 
ential cross sections calculated in the rest frame of the 
emitting source to the laboratory system when compar- 
ing the model predictions to experimental data. It can 
be shown that the transformation may be performed by 
following formula: 



dEdVL 



P_ 

P 



dE'dn' 



d 2 a 

dE'dn' 



where the first equality is exact and the second is valid 
in nonrclativistic limit, normally realized in the motion 
of observed ejectiles. The nonrclativistic relationship be- 
tween kinetic energy E of the particle emitted at the 
angle 9 lab in the laboratory system and the energy E 
of emitted particle in the rest frame of the source is as 
follows: 



E = E 



m/3 2 



\/2mE ■ ■ cosOlab 



where m (= Ap) is the mass of the emitted particle and 
(3 - the velocity of the source in the laboratory system. 
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